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Aspects of meson properties in dense nuclear matter 

Octavian Teodorescu, Abhee K. Dutt-Mazumder* Charles Gale 
McGill University, Montreal, Quebec, Canada 

We investigate the modification of meson spectral densities in dense nuclear matter at zero 
temperature. These effects are studied in a fully relativistic mean field model which goes beyond 
the linear density approximation and also includes baryon resonances. In particular, the role of 
A*(1520) and A*(1720) on the p meson spectral density is highlighted. Even though the nucleon- 
nucleon loop and the nucleon-resonance loop contribute with the opposite sign, an overall reduction 
of p meson mass is still observed at high density. Importantly, it is shown that the resonances 
cause substantial broadening of the p meson spectral density in matter and also induces non-trivial 
momentum dependence. The spectral density of the ao meson is also shown. We study the dispersion 
relations and collective oscillations induced by the p meson propagation in nuclear matter together 
with the influence of the mixing of p with the ao meson. The relevant expression for the plasma 
frequency is also recovered analytically in the appropriate limit. 

PACS numbers: 25.75.-q, 25.75Dw, 24.10Cn 


I. INTRODUCTION 

Electromagnetic radiation constitutes a privileged probe of the properties of matter under extreme conditions. This 
owns partly to the fact that it decouples from the strongly interacting system without significant rescattering and also 
because the virtual photons enjoy a direct coupling to vector mesons. Lepton pairs carry thus valuable information 
about the in-medium properties. 

Among the light vector mesons, the p acquires a special importance because of its large decay width. Therefore, 
this might serve as a chronometer and thermometer to probe temporarily produced hot and dense hadronic matter. 
Even though w or </) mesons do not have this desired short lifetime, in the medium they might undergo sufficient 
broadening leading to a shorter lifetime Q. For the present purpose, however, we first mostly concentrate on the p 
meson and we discuss the scalar-isovector sector later on. 

The in-medium properties of the p meson have been estimated in variety of models ranging from QCD sum-rules 
1^ to chiral models like Nambu-Jona Lasinio, effective hadronic Lagrangian approaches [|j and mean-field models. 
It is fair to say that, at this point, a clear consensus is still lacking but that important progress has been realized 
over the past few years. For a review see 0. The angle of this work consists of uniting several physical aspects 
we find important but that had not been treated together in a unique approach. So here, the p spectral density in 
dense nuclear matter is studied and the importance of the A^*(1520) and Al*(1720) is reiterated, but a relativistic 
calculation is performed and we go beyond the linear density approximation (LDA). A relativistic calculation has 
recently been presented in Q in the LDA. We show that LDA is a good approximation for densities below nuclear 
matter density, but for higher densities multiple scattering becomes important. We report a quantitative comparison 
of the results obtained in the linear approximation with a resummed one loop calculation. We also incorporate 
the effect of interacting nuclear matter through the scalar and vector meson mean fields, motivated by the Walecka 
model. The role of resonances and that of the nucleon loops are examined separately. Finally, we include the recently 
discussed mixing effects and report on the spectral density of the oq for the first time. 

The paper is organized as follows. First the formalism is outlined followed by a discussion of the p meson properties 
involving nucleons. Then we consider the effect of the resonances on the in-medium spectral densities. Later we discuss 
the effect of mixing. We also present the spectral density of the ao meson which supplements our understanding of 
the mixed propagator of the p in nuclear matter. The calculations are done in a fully relativistic formalism including 
the effect of the mean field. At places the mathematical details are relegated to the appendix. Finally we discuss the 
results and conclude. 
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II. FORMALISM 


It is well known that the spectral density is actually the imaginary part of the propagator which in turn is related 
with the polarization functions of the meson. Therefore, we first discuss the properties of the p meson polarization 
function in dense nuclear matter. 

Essentially, the spectral density is related with the collective excitation induced by the p meson by its propagation 
in nuclear matter. This is analogous to the photon propagation in QED plasma where the propagating particle picks 
up the collective modes from the system arising out of the density fluctuation. This is commonly known as plasma 
oscillation. Even though in the present work the main focus is not to recover the characteristic features of the plasma 
oscillation induced by the p meson, nevertheless, we outline the formalism a la Chin ||^ in order to be able to discuss 
the spectral density in terms of the dielectric response function of the nuclear matter. This enables us to incorporate 
the effect of meson mixing in a straightforward manner ^. 

The p meson, being a massive spin one particle, can have both longitudinal and transverse excitations depending 
upon whether its momenta is perpendicular or parallel to the spin. Furthermore, in matter the Lorentz symmetry 
is broken and these two modes will have different characteristic features. These states are designated as n^((7o, |q|), 
n^(9oj |q|) with L and T denoting longitudinal and transverse modes. 

As already mentioned, we consider the coupling of p meson with n-n, n-R and tt-tt states and therefore what we 
have is the following 


= + (i) 

with R = A*(1520), A^*(1720). First we present a general formalism without the effect of mixing and later we shall 
address the issue of the possible mixing and of the corresponding modifications. 

To describe the nuclear matter ground state we invoke the mean field approach of quantum hadrodynamics and 
consequently the effective nucleon mass is generated through the a meson tadpole of scalar mean field potential . 
The nucleon mass is determined by solving the following equation self-consistently. 




d^k ml 

rria Jo (27r)3 ^*2 ’ 


( 2 ) 


To study the collective excitation of the system, the relevant quantity is the dielectric function which actually 
characterizes the eigenvalue condition for the collective modes. In the language of field theory this is equivalent to 
solving the Schwinger-Dyson equation to determine the dressed propagator. The relevance of summing over the ring 
diagrams for the study of vector meson propagation is discussed at length in Ref. |§]. 

The vector meson propagation is calculated by summing over ring diagrams, a diagrammatic equivalent of the ran¬ 
dom phase approximation (RPA), which consists of repeated insertions of the lowest order polarization, as illustrated 
in Fig.^ [|| below. 



FIG. 1. Ring diagrams relevant for the random phase approximation. 

We make use of the Dyson equation to carry out the summation 

= KM + (3) 

The poles are found from the equation 

detys'; - = O. (4) 

The bracketed term is nothing but the dielectric tensor of the system 
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< = ( 5 ) 

the determinant of which, denoted later by e{q), is the dielectric function. The eigenconditions for collective modes 
can now be expressed as e{q) = 0. The relevance of the set of ring diagrams and the origin of such an eigencondition 
can be understood from linear response theory where the fluctuation of the current density, the source term for the 
meson field in nuclear matter, is “picked up” by the vector field. 

For later convenience we define longitudinal and transverse dielectric functions as 

CLiq) = (1 + Z?°noo)(l - D°U33) + 

eT{q) = 1 - D°nii = 1 - D°U22 = 1 - D°Ut. 


( 6 ) 

( 7 ) 


Here, the modified H functions are used and the Dirac part has also been incorporated. Dq = l/{q'^ — m^) is the free 
vector meson propagator of mass my. The eigenmodes of the collective oscillations are given by 

e{q) = e%{q)eL{q) = 0, (8) 

corresponding to the degrees of freedom of a massive vector particle. The two identical (or degenerate) transverse 
collective modes are each given by 


and the single longitudinal mode by 


eriq) = 0 


eiC?) = 0, 


(9) 

( 10 ) 


which yield the relevant dispersion curves. 


A. Pion-pion loop 

It is well-known that the free space decay width of the p meson is dominantly determined by the two-pion channel. 
In other words its coupling to tt — tt loops determines the shape of the free space p spectral density. 

The interaction between a neutral vector meson and the pions are given by 

= gpT.T.i'K X dpTi) ■ p^. (11) 


K 
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FIG. 2. p — rnr Loop. 


The real and imaginary part of the pion-pion loop have been discussed at length in many places |11|J12(| , we just 
quote the results here: 


^gn^(T) _ ^ 4 to^ 3/2 

" 48^2 ^ 2 ) 


1 + t/l- 


1—4m2 
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1- \/i- 


1—4m2 
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( 12 ) 


= - 


/OTTTT 


m2 




487r ' M2 

The free-space spectral density of p meson is given by the following expression: 


Sp{q^) 


_ ImY,p{q^) _ 

TT (g2 — _ Y,pY + ImYij, 


(13) 


(14) 


In vacuum, the above is a Lorentz invariant quantity and a function of q^. In matter, however, this breaks down and 
we shall have non-degenerate spectral densities for the longitudinal and transverse mode of the p meson. 
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B. Nucleon-nucleon loop 


The p-nucleon interaction Lagrangian may be written as 

An* = 9p[Ni^tN + (15) 

n 



FIG. 3. p — nn Loop. 


n“^ = -tg^vSi J ^Tr[*r“*G(fc + g)*f^*G(fc)] 


where Sj is an isospin factor (S'/ = 2 for symmetric nuclear matter). The vertex for p-nn is 


— 7/j 


2m, 


-cTpi^q 


In Eq. (I^), G(fc) is the in-medium nucleon propagator given by [|13| 

GikoAk\)=GFik)+GD{ko,k) 

with 

+ Wn) 


Gpik) = 


and 


(16) 


(17) 


(18) 

(19) 


Goiko, |fc|) = (^ + ml)—6iko - Eu)e{kF - |k|). (20) 

The first term in G°(fc), namely G^p{k), is the same as the free propagator of a spin i particle while the second part, 

GUk), involving 6{kF — |fc|), arises from Pauli blocking and describes the modifications brought about in nuclear 
matter at zero temperature. It deletes the on mass-shell propagation of the nucleon in nuclear matter with momenta 
below the Fermi momentum. It should be noted that here, unlike the situation obtained in usual held theory in the 
absence of extraneous matter, the standard destruction operator operating on the ground state does not give zero 
as long as the momentum lies below the Fermi momenta simply because there are particles with such momenta in 
the ground state. On the other hand, particles can only be created with a momentum higher than fc/? (because of 
Pauli blocking). This clearly is the origin of the second term in the propagator which has a non-covariant form. 
Actually, the in-medium nucleon propagator is dehned in a preferred frame, the rest frame for nuclear matter, and 
this constitutes another important difference with usual field theory. The Fermi momenta which plays such a central 
role is of course related to the density of baryons through pB = (2fc|.)/(37r^). 

When calculating the polarization function ( p^ ) with the nucleon propagator (p^, there will be terms containing 
“G/’G/’”, “GfG_d + GdGf" and “G/jG//”. The first term accounts for the free part i.e. the contributions of the 
Dirac vacuum (H^,^), while the rest provides the density dependent part of the polarization (H^,^), and we can write 


^pu{q) = 

(<z) 

(21) 


d^kTT[r^GF{k + q)f,GF{k)] 

(22) 

n° (9) - nGsls, J 

d^kTvlr^GFik + q)f,GDik) 

(23) 

+ r^G/ 5 (fc -|- q)T^GF^k) + r^G/ 5 (/c -I- q)T^GF)ik)]- 

(24) 



(25) 
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The free part of the p self-energy denoted by is divergent and therefore needs to be regularized. We used the 


dimensional regularization scheme with the following condition: 

i9"n^(9^)/5(g2)"|M;^M,g2=m2 = 0(n = 0, 1 , 2 ..., oo). 
The real part of the density dependent piece of the polarization is given by 


(26) 


= 


pIttSi 


d^k 


(27r)4 J E*(k) 

T^v{k - q,k) 


5{k° - E*{k))e{kF- I k I) 


Tfj,,j{k,k + q) 


.{k — qY — {k + qY — 


(27) 


Beside this there is another term which involves two 9{kF — |A:|) function. That becomes operative beyond 2 times 
the Fermi energy |^. 

The trace involved in the calculation of the loop diagram has three parts corresponding to vector-vector, vector- 
tensor and tensor-tensor terms in the vertex function ©■ These can be cast into the following form 

Tf. = T;:{k, k + q) + T™(fc, k + q) + Tll{k, k + q) (28) 


k + q) = 4[kf,{k + q)^ + {k + q)f,k^ - k ■ {k + q)gF>^ + 


(29) 


rpVt-\-tV 


{Yk + q)=4.M*^q^Q^, 


(30) 


Klik, k + q) = 16(^)2[Q^,(2(fc • qY - q^k^ + q\k ■ q) - q^- 2q^K^Y- 
Hence the self energy can be written as 

KM = KM + KM id) + KM 


The n^j^(g) functions in this case are as follows 


n 


vv 


j-rvt-Ytv 


n 


tt 

flU 


-QM■ pf 

7r3 Ko E*{k) q‘^-4{k-qY 



d^k 1 
E*(k) q'^ — 4(A: • qY 
d^k 

E*{k) q'^ — Ai^k ■ qY 


(31) 

(32) 

(33) 

(34) 

(35) 


where = {k^ - ^q^lik^ - ^qv), Q^i' = (-Pf’' + E^ = Vk^ -|- M*2. To include the overall degeneracy 

factor, the above expressions are multiplied by a factor of two coming from the nucleon and proton loop. It is clear 
that the form for the polarization tensor conforms to the requirement of current conservation, i.e. 

q^K^ = 0 = KM- ( 36 ) 

In the present case we observe that the free part and the dense part of the polarization tensor individually satishes 
the above condition. 

We should also observe that Eq. (|4|) is proportional to and therefore contribute equally to the longitudinal 
and transverse mode. In fact, it is K^i, which in matter, induces the splitting of these two modes as we shall discuss 
later. Evidently, the Dirac part is also proportional to and therefore the modes remain degenerate on account 
of Lorentz symmetry. It might be mentioned that at |g| =0 they are degenerate in matter because of the rotational 
symmetry. 
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Also, it is worth to point out that we could describe these effects in the linear density approximation for low baryonic 
densities. In this approximation, the term cannot contribute and therefore the longitudinal and transverse part 
become degenerate and Eqs. (^3[]^) give 


HtTl = -^9l 


M*ql 




PB 


n 


vt+tv 

T{L) 


2,kM* 


yl A J ) 


n: 


T(L) 


= 4qA_ 


AM ' - AM*^q^) 

M*q‘^ 


PB 


2M' q*-AM*^ql 


PB- 


(37) 

(38) 

(39) 

(40) 


These results could also directly be obtained by multiplying the forward scattering amplitude with the density. 
In the present case we implemented the same by expanding the integrand up to second order or in other words 
schematically we retained terms like d^kf^i,{qo, |(;|)(1 + 0{k)). 

To provide further insight, one can make a long wavelength approximation and recover the familiar results of the 
plasma oscillations. In this limit, when go < Ep and |g| < kp, Eq. ( |^ ) reduces to = gl/M*pB- In this limit, 

with K = 0 the dispersion relation of the density dependent part alone becomes 


ql = \q\^ + ml + n^, 


(41) 


where the plasma frequency = gy/M*pB- This is the non-relativistic form of the results presented by Chin ||] 
for the case of uj meson propagation in nuclear matter. Furthermore, replacing g^ by the electronic charge “e” and 
putting rriv = 0 what one obtains is nothing but the familiar plasma frequency encountered in condensed matter 
physics ||l4|. 

In Fig. ( 3 ) we compare resummed one loop results of the p-meson self-energy with the ones calculated in the linear 
density approximation. The ratio of the self-energies are shown as a function of density. 



FIG. 4. Comparison of the p-meson self-energy calculated at the one-loop level with the one calculated in the linear density 
approximation. 


It is now apparent that the linear approximation can be justified only up to normal nuclear matter density and we 
notice a rather strong dependence on energy, the approximation becoming worse for high values. 
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C. Nucleon-Resonance Loop 


In the present work, among the resonances, we consider only iV*(1520) and A^*(1720) which couples strongly with 
the p meson as indicated in |15| . The corresponding relativistic interactions are given by 

f for /(J-) = i(f) 

Ant = < ' (42) 

for /(J-) = 1(3 + ) 

Here t/ji* denotes the resonance spinor and ip the nucleon spinor, | [ 7 +, 7^] and p'^ — p^. It 

might be mentioned that Ai*(1520) and iV*(1720) have total widths of ~ 120 MeV and 150 MeV respectively with 
corresponding branching ratios into the p-N channel of ^ 24 MeV and 115 MeV . 


k k 




FIG. 5. p self-energy. 


From the Fig.^ it is evident that the polarization tensor has two parts designated as direct and exchange term. We 
present only the direct term and the exchange term can be calculated in a similar fashion. 


- = Si 


d'^k T^^{k,k-q) 


[1 


+ - ml)0{kF - |fc|)], 


(27r)'l (fc — g)2 — 77r|j Lfc2 _ ^2 ' 

where Sj is the isospin factor. The vertex factors for R^^'^Np could be written as 

( 75 ) ^(7/.9a - !^9A*a) 


(43) 


(44) 


It is evident that like the N-N loop, also contains a “free” and a density dependent part. The detailed 

expression for the free part is given in the appendix and has a form i,n(q^). Note also that 

Ti^Ak + q,k)= Tr[i{^ + 


= Tr[{^ =F TO„)(7^(j'a - i9iLa){^ -i + mu) 


- q){q0iu - igpv)]- 

(45) 

In the above equation TVAiAkP) is the Rarita-Schwinger propagator, given by 



(46) 

, 1 M lp+7"-7>^ 

= {p^mR)[ + 7^7 + . 2 q 

0 0 3 rriR 

(47) 

There is an overall sign ambiguity with spin 3 A particles which arises from the special choice of the point transfor¬ 
mation properties of the spin 3/2 Lagrangian |l 6 |,|l^. This has also been discussed in Ref. [^. However, we do not 
adopt here the prescription suggested in Ref. 0 to deal with the spin 3/2 propagator. 

The relevant trace can be written in the following suggestive form 

T^Ak, k-q) = ( A) 2 „±q^^ + 

TTlp 

(48) 


where. 
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( 49 ) 


„±_°r l 2„2 , _ __ 2 , k'^q^ _^k‘^q^{k-q) q^{k-q) , , k'^{k-qf , ^q^k ■ {k ■ qf 


a = - [-k q + runmiiq + 


= i^{k-q-k^- m^)-^. 


ik-qY 


^2 21 
+ 2 - 


•-R 


^R 


It is evident that this structure is similar to what we had for the nucleon loop and therefore satisfies the condition 
of current conservation = 0 = Hf^^q'^ in the momentum space). involves the same gauge invariant forms 

^flu and Qflu ■ 

In order to evaluate II^j^ conveniently, we choose q to be along the z axis i.e. q = (go, 0,0, | q |), and k ■ q = 

E* {k)qo ~ |fc| \q\Xi where y is the cosine of the angle between k and q. After (j) integration the non-vanishing components 
are 


/ IIoo 0 0 Iloa 

0 Hii 0 0 

0 0 1122 0 

V IIso 0 0 Has 

Moreover for isotropic nuclear matter we have Iff^ = and 11^]^ = nf(), and hence, taking all this into account, we 
have only two non-vanishing independent components of linear combinations of which gives us the longitudinal 
and transverse components of 11^,,, in other words: n£(g) = —11^ -|- n|^ and n.^(g) = = 11^. 

We can now estimate the dispersion curves of the p meson in nuclear matter. As mentioned, they appear as poles 
in the propagator and therefore zeros of the dielectric functions shown in Eqs. (^ and In Fig. ^ we show the 

dispersion curves for p meson for p=1.5po baryonic density. It is clear that the p meson physical mass (when qz=0) 
drops in nuclear matter from its free space value. The dashed curve shows the results with n-n loop only and the 
solid one corresponds to the case where the direct coupling of the p with A^*(I520) and Af*(1720) is also considered. 
It should be noted that the resonance-particle excitations contribute with opposite sign to that of the n-n loop. This 
partially offsets the dropping of p meson mass in nuclear matter. 




FIG. 6. The dispersion curve for p showing explicitly the effect of the baryonic resonances. See the text for details. 

Fig. 0 shows the variation of the invariant mass of the p meson mass as a function of nuclear density. The Dirac 
vacuum and the density dependent part of the self-energy contribute with opposite sign to the invariant mass. 


8 













FIG. 7. The invariant mass of the p meson as a function of density showing explicitly the effect of the baryonic resonances. 

At lower densities the free part is responsible for the dropping of the p mass while at higher densities the masses 
again tend to increase when the contribution from the density dependent part increases. This behaviour was also 
observed in the case of the to and a mesons in Ref. [^. 


D. p-ao Mixing via N-N Loop 

Before we start the discussion on the p-ao mixing involving n-n polarization in nuclear matter, we should say a few 
words on the oq coupling to the nucleon. A more detailed study on the oq propagation in a dense medium could be 
found in . The interaction can be described by the following Lagrangian: 

^int — ,a^ 

where ^ and (pao correspond to the nucleon and oq fields, and Tq is a Pauli matrix. The values used for the coupling 
parameters are obtained from Ref. 0 - We do not involve the coupling of oo to the baryonic resonances since currently 
this is not precisely known. 

The polarization vector through which the oq couples to p via the n-n loop is obtained by evaluating the Feynmann 
diagram depicted in Fig. ^ and is given by 


n 


^0 



n 


FIG. 8. ao-p mixing via nucleon-nucleon loop. 


n;.(<3'0, M) = 2igao5p 




Tr[G(fc)r^G(fc + g)], 


where 2 is an isospin factor. With the evaluation of the trace and after a little algebra Eq. 
suggestive form: 


(52) 

can be put into a 
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n^(go, kl) 


9p9ao 


2q\2ml 


2mn Jo E* {k) q^ - A{k ■ qY 


(53) 


This immediately leads to two conclusions. First, it respects the current conservation condition, = 0 = 

Secondly, there are only two components which survive the integration over azimuthal angle. This guarantees that 
it is only the longitudinal component of the p meson which couples to the scalar meson while the transverse mode 
remains unaltered. Furthermore, current conservation implies that out of the two non-zero components of 11^, only 
one is independent. 

In presence of mixing the combined meson propagator can be written in a matrix form where the dressed propagator 
would no longer be a diagonal matrix: 


V = V° + V°IW. 

It is to be noted that the free propagator is diagonal having the following form : 

0 


= 


0 Ao 


In Eg. (|55|) the noninteracting propagator for oq and p are given respectively by 


^o{q) 

D%{q) 


1 


9^ - ’ 



q'^ — rin?p + ie ’ 


(54) 


(55) 


(56) 

(57) 


In fact, it is the polarization matrix which involves non-diagonal elements as shown below characterizing the mixing 


n = 


( n^.(9) 
I nM(9) 


n.(g) \ 

n-o{q)J 


After (j) integration the non-vanishing components 11 are as shown below 


o 

o 

0 

0 

Hoa 

Ho \ 

0 

nil 

0 

0 

0 

0 

0 

n22 

0 

0 

IIso 

0 

0 

Has 

Ha 

V Ho 

0 

0 

Ha 

n“o/ 


For oo meson the free part of the self-energy is given by : 

3o2 

n““(g^) = °° [3 (to*^ — nY) — 4(m* — m)m — {m*^ — rtY) 


— dx {vY — x(I — 

Jo 


dx In 


x)q‘^) In 


— x{l — x)q^ 
rY — x{l — x)q^ 
— x{l — x)q^ 
rY — x{l — x)q^ 


To determine the collective modes, one defines the dielectric function in presence of the mixed terms 


e{qo, Idl) = det(l - I?°n) = x ernix 


(58) 


(59) 


(60) 


(61) 


where er corresponds to two identical transverse (T) modes and Cmix correspond to the longitudinal mode with the 
mixing. The latter also characterizes the mode relevant for the oq propagation 

ct = 1 — doIlT, do = ^ 

r - 

Y 

^mix = (1 ~ doIlL))! — AqIIs) — Aodo(no)^. (62) 

W 

It is evident that only the longitudinal component gets modified because of the mixing, and when Hq = 0 we recover 
the same expression as Eg. dl). 
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III. p SPECTRAL DENSITY IN NUCLEAR MATTER 


Unlike in free space, the longitudinal and transverse p spectral densities are non-degenerate. This is because of 
the fact that, in the nuclear matter rest frame, Lorentz symmetry is broken and they are functions of go and |q|, 
independently. Furthermore, in matter the scalar and vector mesons can mix. This also modifies the longitudinal p 
spectral density through the off-diagonal mixing terms in Eq. (^). 

Now in the presence of mixing, the spectral densities can be defined in terms of the dielectric function as 

SL{qoM\,PB) = --Im[^^^^^^^]. (63) 

TT esL 

On the other hand the transverse spectral density is unaffected by the mixing and has the following form 

5'T(go, |q|,PB) = (64) 

TT i — CLqLLx 




M. (GeV) 


FIG. 9. The longitudinal spectral density for the p meson with mixing at density p = 1.5po- 

Fig.| shows the longitudinal spectral density (Sl) at a density p = 1.5po as a function of three momenta (|q|=g 2 ) 
and invariant mass (M^). This includes the effect of N-N loop and the direct coupling of the p meson with iV*(1520) 
and A* (1720). The three-peak structure of the spectral density is similar to what has already been observed in non- 
relativistic calculations in Ref. [Q . Such a characteristic behaviour of the spectral density in presence of resonance 
has been also observed in pion-nucleon dynamics. The collective modes induced by the density fluctuations can be 
identified as the p meson mode, N-hole and N-resonance mode similar to what we observe in case of pion propagation 
in nuclear matter. 
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FIG. 10. The longitudinal spectral density for the p meson with mixing at density p = 2.5po. 

Fig. 1^ is as Fig. ^ but at higher density. We note that at higher density the spectral density gets even broader. 
Furthermore, with increasing momenta all the peaks merge into one broad peak indicating the fact that at high 
momenta the collective behaviour dies down and the meson shows a behaviour of a free propagation in matter. 



FIG. 11. The longitudinal spectral density for the p meson with mixing at density p = 2.5po with and without considering 
resonances. 

Next we show how individual resonant states modify the p spectral density in matter. In Fig. |^, the dashed-dotted 
line represent the p spectral density for p coupled only to tt-tt and N-N loop. This shows the p meson peak shifts 
towards lower invariant mass. 
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FIG. 12. The transverse spectral density for the p meson with mixing at density p = 1.5po. 

This is because of the a meson mean field in nuclear matter which lowers the nucleon mass considerably in nuclear 
matter. The dashed curve has a double hump structure with the introduction of the Ai*(1720). For non-relativistic 
calculation, such a feature has already been observed in the static limit in Ref. |]^. Addition of Ai*(1520) also 
introduces a double peak structure but it moves the strength to the low invariant mass region. Hence we conclude 
that the resonant states are responsible for the broadening of the p spectral density in nuclear matter. Furthermore, 
in presence of the resonant states the p meson gets broadened so much so that the single particle interpretation of p 
meson as quasi particle fails to carry any sense. 

The transverse spectral density {St) also shows interesting features at low momenta. Again we observe that at 
higher momenta the p-spectral density gets flattened. The momentum dependence of St is also observed to be 
different than that of ^l. 



1.2 


FIG. 13. The transverse spectral density for the p meson with mixing at density p = 2.5po. 
At higher densities we see even more pronounced effect as evident from Fig. |^. 
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The spectral density for the oq meson is defined as: 


Of II ^ ^ T r^o(l — dolls)^ 

Ss{qo, |q|, Pb) = - Im[ -] 

71" esL 


( 65 ) 



FIG. 14. Spectral density for the ao meson with mixing at density p = 1.5po. 

In Figs. 0 and 1^ the spectral densities of the oq meson are presented at densities p = 1.5po and p = 2.5po- It is 
evident that the Op spectral density also gets broadened in medium. A marked shift towards the lower invariant mass 
indicates that the op mass also drops in nuclear matter. 



FIG. 15. Spectral density for the ap meson with mixing at density p = 2.5po. 

To illustrate the spectral density modification in matter, we present 5^ as a function of density. The solid line in 
Fig.|l^ shows the free p spectral density. We find that with increasing density it acquires more strength in the low 
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invariant mass region and also becomes flattened. Dashed, dotted and dashed-dotted lines represent results for 0.5, 
1.5 and 2.5 normal nuclear matter densities. It might be mentioned that a small peak appears at the higher invariant 
mass region which indicates the effect of mixing. 




FIG. 16. The longitudinal spectral density for the p meson with mixing as a function of density. 
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FIG. 17. The spectral density for the ao meson with mixing as function of density (upper part) and a cross section for 
different densities (lower part). 

In Figjl^, the two dimensional projection of the oq spectral density at |q| = 0.7 shows that with increasing density 
ao gets narrower and moves towards low invariant mass region. This is understandable from the reduction of phase 
space. It might be recalled here that resonant states do not couple to ao and therefore we do not observe any 
broadening in the oq channel. 


IV. SUMMARY AND CONCLUSIONS 

In the present work we present quantitative results of the in-medium meson properties. The meson spectral densities 
are evaluated for the first time in a fully relativistic mean field model which goes beyond linear density approximation 
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and we discuss the limitation of the aplicability of LDA. We find that meson spectral densities in matter are quite 
different than in free space. The difference stems partly from the Lorentz symmetry breaking. While in free space 
the transverse and longitudinal component of the p spectral density are degenerate, in matter, they show different 
qualitative behaviour. Furthermore, we include other purely in-medium effects, forbidden in vacuum on the account 
of Lorentz symmetry, like meson mixing. In matter the p can also be modified because of its mixing with a scalar 
(isovector) oq meson. This effect, as we have seen, modifies only the longitudinal component of the spectral density. 

In our model we observe that the p meson spectral density gets flattened in nuclear matter with the incorporation 
of the resonant states like 7V*(1520) and A^*(1720). In fact in nuclear matter the original free space peak of the p 
meson become so broadened that it is no longer possible to interpret that as a quasi-particle excitation. This was also 
observed in Ref. [pAplj . It should also be noted that in presence of scalar mean field the p meson mass goes down as 
a function of density. 

This broadened of p spectral density shows an accumulation of strength towards lower invariant mass region. This 
would definitely imply that in more production of dilepton pairs with low invariant mass in access to what we might 
expect from the free p spectral density. We plan to extend this work to the finite temperature region in the near 
future. Furthermore, the broadened p spectral density would also affect the width of the resonance states like A^*(1520) 
or A^*(1720) which probably would shed some light on the issue of mixing resonant states in the photoabsorption 
cross-sections. For the complete determination of the p spectral density and the broadening of the resonant state, a 
self-consistent approach should be adopted Q. Studies along such directions are in progress. 
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APPENDIX A: FREE PART OF THE N-N LOOP 


Polarization tensors arising out of the n — n excitation of the Dirac sea. 


n 


vv 

fiiy 


IT o 



— x)lnD\ 


(Al) 


xrRi+iR 


1 gl M*ti, 

2 7r2 2M 


Qf^^[A - 1 - Inp,'^ 



(A2) 


Rtt 


(47r)2'M 


rQ, 


[ dxx^lnD — q^ [ dxxlnD — [ dxlnD) (A3) 

6 o2 ./n ./n ./n 



7 -|- InAi: 


D = -q^x{l-x) 

All the terms containing A are infinite and need to be subtracted out. 


APPENDIX B: FREE PART OF THE N-R LOOP 

J^RN{d^r) ^ Q^^n(g2) 


where. 


^(g^Ai) = (ni(g^) -Hn2(g^) -f n3(9^) -tn4(g2) + n5(g^) -tn6(g^) -f n7(9^) -tn8(g^)) 


(Bl) 


(B2) 
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(B3) 


Ui{q^) = [{D + 2x^q^){A + \nm^)+ [ dx{D - D\nD - 2q‘^x‘^\nD)]'^^;^ 

Jo 327r 

U2{q‘^) = [2(3i:)2 + 18i:)gV+4gV)(A + ln/i2) + 

+ / dx{9D‘^ + 36Dq^x^ + 6D^ InD — 36Dq^x^ \nD — Sg^a:^ In D)] ^ ^ 

Jq IZoTTlp^Tr 

113 ( 9 ^) = -[2g^a;(3L) + 2g^a:^)(A + In^^) + 

+ f dx{2>q^xD — Sq^xD InD — 2q‘^x^ In D)] ^ ^^‘'2 

Jq OZTTlj^Tr 

114(9^) = -[{ 2 D + q^x^){A + In/i^) + 

+ [ dx{D-2D\nD-q^x‘^\nD)]%^ 

Jo IStt 

n 5 (' 7 ^) = -TOArmA[(A + In/i^) - y dx In D] ^^^2 
Ue{q^) = [D{A + In^i^) + j\x{D - DinD)]^^ 

Ur{q^) = [2D{3D + 2q^x'^){A + \nfi‘^)+ f ^D{9D + 4q'^x'^ - 6D\nD - Aq^x^lnD)]'^^^ 

Jo TOjj 1287r"^ 

n8(g^) = -[ [ dxDq‘^x{A + lnf/) + [ -^{Dq^x - q"^ x D In 

Jo Jo oZTr 


APPENDIX C: IDENTITIES INVOLVING RARITA-SCHWINGER SPINORS 


1 2 

^tf + M«)[-9- + 37'‘7" + 5^ 
= [t + M„,)P^J,(p) 


1 


3 Mr 


(Cl) 


Clearly the above and the following equations vanish when the spin 3/2 state is on-shell, so does Eq. (C3) also. 

Another useful relation in this respect is the identity 

2 (r)^ — 1 2 

{k-p)^P^j2{p){k-p)^ = - ^2 ^'^ [p^-2fc-p] -h ■ pf - Mlk"^] (C4) 

It is to be noted that the above equation also takes much simpler form when we have on-shell spin 3/2 projection 
operator. 

p^^'^(p)=0 = 4'"(p)p, (C5) 

k^P^J^{p)K = ^[(^ • Pf - Mlk^] (C6) 

It is evident that Pfj.Py 2 {p)Pi> = 0 when the spin 3/2 particle is on-shell i.e. p^ = M|.. In the rest frame of the 
resonant state R, we, therefore, have kfj^Pl^J'^kv = |k^- 
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